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Abstract 

The conventional theory of superconducting alloys does not take 
into account a discreet character of impurities. Experimental data for 
superfluid 3 He in aerogel and for some of high- T c superconductors 
reveal a significant discrepancy between the observed temperatures 
T c of their transitions in the superfluid or superconducting state and 
that predicted by the theory. Here a theoretical scheme is presented 
for finding corrections to the T c originating from spatial correlations 
between impurities. Analysis is limited to the Ginzburg and Landau 
temperature region. The shift of T c with respect to the pure material 
is represented as a series in concentration of the impurities x . In the 
first order on x the conventional mean-field result for the lowering 
of T c is recovered. Contribution of correlations enters the second or¬ 
der term. It is expressed via the structure factor of the ensemble of 
impurities. For superfluid 3 He in silica aerogel the sign of the correc¬ 
tion corresponds to an enhancement of the T c so that the resulting 
pairbreaking effect of impurities is weakened. When correlation ra¬ 
dius of impurities R exceeds the coherence length of the superfluid 
£o the contribution of correlations to the shift of T c acquires a factor 
~ (i?/£o ) 2 and the weakening of the pairbreaking effect becomes ap¬ 
preciable. The presented scheme is applied to the superfluid 3 He in 
aerogel. 
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1 Introduction 


Impurities can significantly change properties of a superconductor. E.g. it 
is possible to change a type of superconductor from first to second by an 
increase of concentration of impurities. In case of unconventional Cooper 
pairing impurities lower transition temperature T c . The word ’’unconven¬ 
tional” means that at this type of Cooper pairing not only gauge symmetry 
but some other symmetries of a normal phase are broken. Transition tem¬ 
perature is one of the most important thermodynamic characteristics of a 
superconductor. It is also one of the first quantity to be measured in ex¬ 
periments. Abrikosov and Gorkov (AG) [T] theory of superconducting alloys 
renders quantitative description of the lowering of the transition temperature 
and of other properties of superconducting alloys in terms of one parame¬ 
ter - a transport mean free path l tr of electrons in a given material, which 
in its turn is inversely proportional to the number of impurities per unit 
volume n , or their concentration x . In the AG-theory it is assumed, that 
distribution of impurities is completely random, i.e. their positions are not 
correlated. Effect of correlations on thermodynamic properties of supercon¬ 
ductors is of the second order on concentration and in most cases it can be 
neglected in comparison with the principal mean field effect, which is linear 
in x , but there may be special reasons for an enhancement of contribution 
of correlations. 

In the last two-three decades many unconventional superconductors were 
discovered, including liigh-T c materials. Doping with impurities is signif¬ 
icant part of a processes of their preparation. The pairbreaking effect of 
impurities in high- T c superconductors is observed experimentally M, but 
the lowering of T c in many cases is smaller than that, predicted by the 
AG theory (cf. [I] and references therein). This tendency is interpreted in 
Ref. [4] as a manifestation of inhomogeneity of the condensate which can be 
appreciable in superconductors with a very short coherence length. 

Doping of metallic superconductors changes not only concentration and 
distribution of scattering centers but also other parameters like concentration 
of charge carriers, lattice constant etc., which in their turn change the T c . 
A special effort has to be paid to separate effects of different factors. The 
mentioned difficulties are not present in another physical object, where effect 
of correlation of impurities is observed - it is superfluid 3 He in aerogel. 
Liquid 3 He becomes superfluid at temperatures of the order of 1 niK. 3 He 
is unconventional superfluid - Cooper pairs here are formed in a state with 
angular momentum l = 1 and spin s = 1, so that the rotational symmetries 
both in spin and in the orbital space are broken. Properties of pure (or bulk) 
superfluid 3 He are thoroughly investigated experimentally and most of them 
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are quantitatively interpreted [jS] so that changes introduced by impurities 
can be accurately separated. All that makes liquid 3 He a favorable object 
for a study of effect of impurities, except that introduction of impurities 
in the superfluid 3 He is not a trivial task. Possible floating impurities at 
these low temperatures stick to the walls of a container. To remain in the 
bulk of the liquid impurities have to form a self supporting structure. As 
this in experiments high porosity aerogels are used PEli. Majority of 
data is obtained with silica aerogel. It consists of very thin (diameter d ~ 
3 nm) strands of SiO 2 ■ Scattering of quasi-particles by the silica strands 
renders a finite mean free path l tr and lowers a temperature of transition 
in the superfluid state. Instead of concentration aerogels are traditionally 
characterized by the porosity P , i.e. by the fraction of space which is not 
occupied by the strands. In most of the experiments with 3 He silica aerogels 
with P cs 98% are used. In theoretical models the strands forming aerogel 
are often treted as chains of silica balls of a radius p ~ d/2 , which play the 
part of scattering centers. The centers can be treated as impurities in the 
conventional theory as long as p <C £0 where £0 is the coherence length of 
the superfluid. 

Essentially new property of aerogel is that for forming a rigid structure po¬ 
sitions of the centers have to be correlated in space. The correlation function 
is an intrinsic property of aerogel, it does not depend on a substance filling 
its pores, and it can be directly measured by X -ray scattering. Presence 
and importance of correlations was realized at the very beginning of investi¬ 
gation of 3 He in aerogel [9j. I 11 particular correlations were suggested to be 
responsible for deviations of the observed dependence of on y from that 
predicted by the AG theory. For curing the disagreement more complicated 
theoretical models of 3 He in aerogel were suggested uni- One of them com¬ 
bines description in terms of impurities with that in term of pores (isotropic 
inhomogeneous scattering model IISM). This model improves agreement with 
the data for 3 He but it has certain drawbacks. The model exploits uncon¬ 
trolled approximations and exact physical meaning of fitting parameters is 
not clear. Another successful model (phenomenological IISM) is an heuristic 
interpolation between the limiting situations of impurities and pores. It is 
not clear also how the IISM can be generalized for application to correlated 
impurities in superconductors. 

In a present paper correlations are treated within the perturbative ap¬ 
proach, developed by Larkin and Ovchinnikov for superconductors with a 
scalar order parameter HU- Effect of quenched disorder is introduced in the 
Ginzburg and Landau equations via position dependent coefficients. Ran¬ 
dom deviations of the coefficients from their average values are treated as 
perturbations. This approach is more standard and more universal than the 
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previously used. A preliminary account of application of such approach to 
the superfluid 3 He in silica aerogel was published before [12] . In comparison 
with Ref. H 2 here a more general and more realistic formulation is used. In 
particular the anisotropy of coherence length of 3 He and possible anisotropy 
of scattering cross-sections of impurities are taken into account. Better esti¬ 
mation of a region of applicability of the present approach is given and the 
process of summation of the perturbation series is fully reconsidered. 

2 Phenomenology 

Pair-breaking impurities make the condensate of Cooper pairs spatially 
nonuniform. In a vicinity of the transition temperature T c free energy of 
such superconductor can be written as a Ginzburg and Landau functional 
with the coefficients, which are random functions of coordinate. For a con¬ 
densate with the scalar order parameter T (r): 

F s {^> (r)} = F n + y'{a(r)|'F(r )| 2 + ife(r)|'L(r )| 4 + c(r)|V'F(r)| 2 }d 3 r. (1) 

In view of further applications magnetic field is not taken into account in this 
functional. For a situation when fluctuations of coefficients a(r), b(r), c(r) at 
their average values are small Larkin and Ovchinnikov m have analyzed the 
effect of inhomogeneities on thermodynamic and electro-magnetic properties 
of a superconductor with the scalar order parameter treating the fluctuations 
as a perturbation. According to their analysis the most strong effect on 
the average value of the order parameter and on the shift of the transition 
temperature have fluctuations of the coefficient a( r), which can be expressed 
in terms of fluctuations of the local transition temperature T c (r): a(r) = 
a(T — T c ( r)). Contribution of these fluctuations is singular at T —> T c . 
Contributions of fluctuations of b( r) and c(r) are regular in this limit, they 
can be neglected so that the average over ensemble of impurities, values ( b) 
and (c) can be used in the equations. If T c0 is the temperature of transition 
in the absence impurities, then a(r) = aT[r}(r) — r], where r = (T c0 — T)/T 
and 77 (r) = (T c0 — T c (t))/T is a random local shift of T c . After introducing 
dimensionless variables 4/ = d'oV’ with d/g = aT(b), and £ 2 = (c)/<yT, or 
£ 2 = '^r^o the equation for the extremum of the functional Eq. (1) takes a 
form 

[T-7?(r)]^ + £ 2 A?/;-V#| 2 = 0. (2) 

Global T c is defined as the highest value of T for which (i(j) 7 ^ 0. The 
signs in the definitions of r and r/ (r) are chosen so that the linear part of 
Eq (1) is analogous to the Schrodinger equation of a particle moving in a 
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random potential r/(r) and r has a meaning of energy. Superconducting 
states correspond to positive r. For small rj( r) a shift of the transition 
temperature r c = (T c0 — T c )/T c can be represented as a perturbation series 
r c = +rj 2 ^ +..., where r^ 1 ' = (?/(r)) and depends on the correlation 

function (r]( r)?/(r’)) . 

The perturbative approach of Larkin and Ovchinnikov can be applied to 
condensates with a multi-component order parameter as well. That requires 
modification of the free energy functional according to the form of the order 
parameter. In superfluid 3 He the order parameter is a complex 3x3 matrix 
\i 0 - Its first (Greek) index corresponds to three possible projections of 
spin and the second (Latin) - to three projections of the orbital momentum of 
a Cooper pair. Non-magnetic impurities interact with the orbital part of 
and the additional term in the density of free energy is , 

where Tjji(r) is a random tensor field. It is assumed that the impurities 
preserve t —> — t symmetry, then rjji(r) is real and symmetric. 

Like in the scalar case, the perturbation can be separated in the ensemble 
averaged part (rjji) = + Kji and the fluctuation fjji( r) = ^(r) — (rjji) 

The isotropic part r 7 ® can be included in the definition of (T c ) so that 
r = an( q j g a global anisotropy. The finite anisotropy splits one 

T c generally into three T c corresponding to different components of l. As a 
result in a neighborhood of T c there may be several transitions into phases 
with different order parameters mm- To avoid irrelevant complications in 
what follows it will be assumed that the ensemble of impurities is isotropic 
i.e. Kji = 0 . The physical meaning of fjji( r) is clear from the definition - its 
isotropic part describes local fluctuations of T c (r), while the anisotropic part 
Tjji(r) = rjji( r) — |-/^(r )5ji describes the local splitting of T c for different 
components of l. With these notations 

Fgl = 

Ew 

S=1 

(3) 

where I s - are the fourth order invariants in the expansion of the free energy 
over A m 0 and j3 s , s = 1, ...5 are phenomenological coefficients. 

For finding of the T c in a presence of perturbation rjji(r) it is sufficient to 
keep in the Ginzburg and Landau equation corresponding to the functional 
(3) only linear with respect to A M terms 

[r Sjl - W (r )]A pl + & + 2 J^) = » (4) 


^(o) / d 3r \ foi(r) - rSj^A^A* + £ 


' dA dl dA*^ 

v dxj dxj 


+ 2 - 


dA jd dA 




dxj dxi 


+ 
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Figure 1: Perturbation series for the Green function of Eq. (4) 

and to use a standard perturbation procedure. Fourier transform of the 
Green function Gji(k, k') of Eq. (4) is developed in the series, graphically 
represented on Fig. 1. Each line corresponds to the Green function of the 
unperturbed Eq. (4) 


G<? ) (r,k) = 


kjk 


T 


3 a# 


+ 




kjki 


T 




(5) 


and each cross - to the Fourier transform of the perturbation rjji (k s+1 — k s ). 
The series has to be averaged term by term over the ensemble of rjji and 
the result is re-summed into the averaged Green function, which is spatially 
uniform: 

(Gji(r, k. k')> = (27r) 3 <5(k - k'JG^T, k). (6) 

In its turn 


G„(r,k) = G‘?(r,k)+Gf2(r,k)(iy m „(k,r))GS(T,k), (7) 

where (W mn (r, k)) is the averaged sum of the series graphically repre¬ 
sented on Fig 1 if the thin lines on both ends of each graph are omitted: 

(Wmn (r, k)) = w%l (r, k) (8) 

p 

If (Wji(r, k)) is small (the criterion will be formulated later), in a principal 
order 

G-/(r, k) = (r - :&)k,h + (r - £* 2 )(fy - kfh) - (W^t, k)). (9) 

Transition temperature r c in a presence of perturbation is found as a pole 
of Gj { (t,0) 

T c Sji = (Wj 7 (t c ,0)). (10) 

In what follows all functions of r will be used at r = r c and the argument 
r will be suppressed. In the first order on the perturbation the mean held 
result is recovered 

^(0) = (hj/( k = 0)) (11) 
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and in the second order we arrive at a straightforward generalization of the 
corresponding result of Ref. m- 

( 2 )/"n\ f ^nm 2fc m A) r j d k f-\r,\ 

W V(°) = -] -- (Vjm(k)r, nl (-k))j^. ( 12 ) 

The variable r does not enter this expression because of the condition (10). 


3 Impurities 


For application of the obtained expressions a form of Tjji(r) has to be speci¬ 
fied. In what follows we assume that perturbation is produced by ensemble 
of discreet identical impurities situated at random positions r a . If linear 
dimensions of each impurity p meet the condition p <C £o they can be con¬ 
sidered as ” small objects” in a sense of the theory of Rainer and Vuorio [15]. 
According to this theory every impurity acts on the condensate as a localized 
perturbation with a characteristic size of the order of £ 0 • For example, if the 
impurity is a ball of the radius p , which diffusely scatters quasi-particles, 
the corresponding perturbation at distance r q from the center is given 
by 



P , . , 
“ v i"‘ 111 



(13) 


where Vj is a unit vector parallel to r. If £ 0 <C l the perturbation pro¬ 
duced by all impurities can be approximated by the sum of contributions of 
individual impurities: 

Vji(r) = ^2v^\r-r s ). (14) 


The Fourier transform of rjji( r) is 


%( k ) = 5>.j?( k ) exp(-?kr s ). (15) 

S 


For substitution in Eq.(ll) we need only Vji\ 0) • For identical 

impurities tensors rff (0) differ only by their orientation. In case of uniaxial 
impurities [15] 




cr 


Sjl + (Ttrfia 


j 



(16) 


where ap is a unit vector in the direction of symmetry axis of the impurity 
at a point r s and the transport cross-sections ap) and cr t “ are expressed 
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through the differential cross-section of scattering of quasi-particles by the 
impurity as: 


- 5 jt ) 



v')- 


(17) 


At the averaging over directions of cip the term, proportional to i n Eq. 
(14) vanishes. Combining Eqns. (8),(9),(11) and (16) we arrive at the first 
order correction to the transition temperature 


Tj? = n ~^£o VfrSji ( 18 ) 

This expression coincides with the result of AG theory in the first order on 
the ratio £ 0 /hr- The correction is negative, i.e. the transition temperature 
decreases. 

To find Tji we have to substitute in the r.h.s. of Eq. (12) the explicit 
form of r/ji(k ). It will be shown that the principal contribution to the integral 
comes from the region of small k . Then the impurities can be considered as 
point-like: 

Vji( r ) = °)<K r - r s) (19) 

S 

with the Fourier transform: 

Vji( k ) = 5] vjf(O) exp(—ikr s ). (20) 


Eventually 



M) 


ti - I d k 38 nm 2k m k n 

~T~J J (27r) 3 3 C s k 2 


8jm.8nl (&tr 


M\ 2 y 


yk(rt—r s )) + ( a (ah 


! (E(3' 


;(«) 


a$ - h jm )(3a«af ) - ^)e ik(rt “ rs) ) 


L M M 

(21) 

If impurities are spherically symmetric =0 ?/jf(0) = rjo5ji. Distribution 
of impurities is characterized by one function - the structure factor 


S( k) = (^e ik(rt - rs) ). 

t 


( 22 ) 


In the structure factor the contribution of correlations can be separated. The 
term with t — s is always present, it is equal to unity and when substituted 
in Eq. (19) renders second order correction to the shift of the transition 



temperature by non-correlated impurities. The remaining sum can be written 
as an integral 

S( k) - 1 = nj C(r t \r s )e-^ rt -^d 3 r t , (23) 

where C(rt|r s ) is the probability to find a particle in the point r t if there is 
a particle in the point r s . For isotropic distribution of impurities it depends 
only on a distance r = |r f — r s | . At r —> oo correlations vanish and C(r) 
tends to a constant. Normalization of C(r) is chosen so that this constant 
is unity. The unity contributes to S(k) a term, proportional to n<5(k) , 
which is already taken into account in r-^ . A measure of correlations is 
v(r) = w(r) — 1. The structure factor can now be represented as S'(k) = 
1 + (27r) 3 n5(k) + S'(k), where only the last term depends on correlations 

5(k) = n J v(r)e~ lkr d 3 r. (24) 


Its contribution to r ^ is: 



5 

2C(3) 



S{k)dk. 


(25) 


The X-ray scattering data for silica aerogels mum show that there exist 
an interval of k: (1 /R) < k < (1/p) where S(k) ~ (1 /k Df ), i.e. these 
aerogels have a fractal structure with the fractal dimension Df ~ 1.7 -t- 1.9 
depending on a sample. The inverse lower boundary of the fractal interval has 
to be identified with the correlation radius R , the inverse upper boundary 
p is of the order of characteristic size of a structure element of aerogel. 
The integral in the r.h.s. of Eq. (25) converges for large k , a principal 
contribution to this integral comes from the region of k ~ (1 / R) . For 
R £o this justifies approximation of rff (k) by 77 ^( 0 ). The available 
experimental data for S(k) are presented in arbitrary units, so that they 

/Q\ 

can not be used directly for evaluation of f>; . Another possibility is to 
express in terms of v(r ), using the relation 

roo _ roo 

/ S(k)dk = 27 t 2 71 / v(r)rdr (26) 

Jo Jo 

and make a plausible guess about the form of the function v(r) , introducing 
in it adjustable parameters. Minimum two parameters are required. One is 
correlation radius R , introduced above, another is the overall amplitude A 
so that v(r) = Av . In this case 




5 7T^ R 2 

C(3) "9" ^ 


(27) 
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where / = / 0 °° v(x)xdx is a number. This number should not be very sensi¬ 
tive to the particular form of v(x) A convenient model expression for v(x) 
in aerogel is mm- 


v(x) 


r (£>/) 


X 


Df—3 



exp(-x), 


( 28 ) 


where T(.D/) is Euler Gamma-function. This expression has a fractal asymp¬ 
totic v(x) ~ x 3 - D at r <C R. At r R it tends to zero and it satisfies nor¬ 
malization condition / v(x)x 2 dx = 0 . For such v(x) I = (3 — Df)/(Df — l) . 
If Df— 1.8 / = 3/2. Experimental data and simulations p] show that in 
aerogel impurities form clusters. This tendency corresponds to Df <3 and 
positive A . In this case the sign of the correction f^ ; is opposite to that of 

Tjl . It means that the destructive effect of impurities is weakened and the 
resulting T c is higher than that given by AG-theory. The increase of T c is 
due to adjustment of the condensate to local inhomogeneities which increases 
a gain of energy. 

The correction fj / 2 ' 1 is obtained with the aid of perturbation theory. It 

means that the inequality |has to be met. This condition sets 
an upper limit for R: ( AR 2 /C, 0 l tr ) <C 1. A limit from below is set up by a 
condition securing dominant contribution of correlations: nAR 2 p 1. So R 
has to be within a window A i. The amplitude A is not measured 

directly in experiments. For its estimation further model assumptions have 
to be used. Assuming, that threads of aerogel are ’’beads” formed by the 
balls of the radius p we arrive at a condition ( 47 r/ 3 )Anh( 3 p ) 3 ~ v , where 
v is a ’’coordination number” i.e. average number of impurities touching 
a selected one. For ’’beads” v « 2. If the fractal dimension Df m 2 
A ~ a(l tr /R) with a ~ (1/10). Within this model we arrive at a stringent 
condition for applicability of a sum as an approximation for r c : 

p -C aR < • Particularly restrictive is the upper limit. Lifting of the 

formulated restriction requires an account of higher order terms in the series 
Eq. ( 8 ). 


4 Long-range correlations 

The series Eq. (8) is an expansion in the parameter Ah- . The p -th term of 
the series 


wjr 1} ( k ) = Vo +1 ™ j 


d 3 ki 

M 3 


d kp 

(2A) 3 jmi 


*(o) 


( k i 


■G®,( k p ) 


10 







( 29 ) 


X) exp[i(k - ki)(r„ - r e )]...exp(t(k p _i - k„)(r 0lJ - r e )] 

ai,...,a p 

Here ai,...,a p ,e are indices numbering impurities. Summation over the 
index e has rendered the factor n in front of the integral. Expression 
in the angular brackets is a product of values of the random function 
5(q) = J2 a exp[?'q(r a — r e )] taken at different values of its argument. Es¬ 
sential contribution to the integral comes from the region k s ~ (1/i?). If R 
is much greater than the average distance between the impurities the num¬ 
ber of impurities contributing coherently to the sum is large. In this case the 
random function S'(q) is close to its average value i.e. to the structure factor 
S( q) = (S'(q)) and in the leading order over nR 3 the result of averaging in 
Eq. (29) can be represented as a product 5(k —ki)S'(ki — k 2 )...S'(k p _i — k p ). 
Then the consecutive members of the series for W (k) are related via 

w^ +1) ( k ) = Vo J - k i)Gjm(ki)wS(ki) (30) 

and summation of the series renders the integral equation for Wji{ k): 

/ (J6U 

^(k - ki)G' jm (ki)H / m/ (ki). (31) 

for substitution in Eq. (10) we need only Wji (k = 0) . The main contribution 
to the integral in the r.h.s. of Eq (31) comes from the region k < (1 /R). 
At W = const, the integral converges, so that Wji( k) fa ^(0)^7 can be 
used as an approximate solution of Eq. (31). For an isotropic W 3 i (k) the 
Green unction Gj m { ki) can be averaged over directions of ki, rendering 
Gjm( k x ) = Then the value of W( 0) is determined by Eq. (31): 



w '-(°> = ivyg. (32) 

where Q = j ^VS'(ki)G(ki). Comparison with Eq.(18) shows, that sum¬ 
mation of the perturbation series reduces to substitution instead of a$ an 
“effective” cross-section, which takes into account correlation of positions 
of scattering centers. For the model correlation function Eq. (28) with an 
account of Eqns. (10),(16) 


jf_to_ _ 1 _ 

4 l tr 1 + 


( 33 ) 


where B(D) 


10r(J>) 3-D 
9C(3) d -i • 
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The obtained expression (33) for a relative lowering of T c has a similar 
structure to that, suggested on a basis of heuristic argument by Sauls and 
Sharma ng. In a limit 1 -Cl r c is proportional to the ’’pairbreaking 
parameter” ff 1 , as in the conventional theory. In the opposite limit the 

c2 

mean free path cancels out and r c ~ R2 is determined by geometric 

characteristics of aerogel. When R is growing r c decreases and T c tends to 
its bulk value. If A is estimated within the ”model of beads” the dependence 
on R is different: r c ~ IRr^bR) d He ^ c anc ^ consequently £0 strongly 
depend on pressure. At solidification pressure £ 0 is approximately four times 
smaller than at the vapor pressure. In 98% silica aerogel both limits of Eq. 
(33) can be reached. With the aid of their interpolation Sauls and Sharma 
P2| were able to fit the pressure dependence of T c in such aerogel for the 
interval of pressures from 6 to 32 bars, using realistic parameters of aerogel. 
Within the region of validity of Ginzburg and Landau equations Eq. (33) 
renders a theoretical justihcation to their interpolation. 

3 He in aerogel was used here as a simple example, where effect of cor¬ 
relations can be analyzed and compared with the existing data. To obtain 
a concrete result in deriving Eq. (33) specific properties of correlations in 
aerogel are used, so it can not be directly applied to a different object. On 
the other hand the argument of previous section, leading to Eq. (31) is gen¬ 
eral and can be used at analyzing of effect of correlated impurities in metallic 
sup er conduct ors. 

5 Discussion 

The example of superfluid 3 He in aerogel clearly demonstrates that corre¬ 
lations of impurities can substantially weaken their pairbreaking effect and 
increase the transition temperature with respect to that expected for non- 
correlated impurities. More literally, it shows that the description of aero¬ 
gel as a uniform continuous media is not sufficient. Correlations introduce 
a spatial scale R. Adjustment of a superfluid condensate to the external 
inhomogeneity increases its gain of energy in comparison with a uniform 
condensate. Mostly important are correlations with a characteristic radius 
R of the order or exceeding the superfluid coherence length £0 • Superfluid 
condensate is more susceptible to such correlations, because reaction of con¬ 
densate on external perturbation is determined by momenta on the order 
of l/f 0 . For comparison, reaction of normal quasi-particles is determined 
by much higher momenta - on the order of Fermi momentum. In this case 
effect of correlations renders only a small correction to the mean free path. 
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Enhanced susceptibility of condensate to long range correlations manifests 
itself below T c as well. A separate analysis shows that in a superfluid phase 
correlations effect a temperature dependence of the order parameter because 
of the interplay of correlation radius and Ginzburg and Landau correlation 
length £(T) fl8l 19]. 

In the last few years experimental activity in 3 He shifted to using different 
types of aerogel. On one hand there are highly anisotropic ’’nematic” aerogels 
0120] on the other - highly isotropic silica aerogels, exhibiting new features 
I ZB- Application of the present approach to these objects can bring further 
interesting results. Practically even more useful and stimulating would be a 
systematic experimental investigation of the effect of correlated impurities in 
metallic unconventional superconductors. 
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